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1. Introduction 

The proper quantification of credit risk poses a complex mix of problems, as important 
credit risk parameters such as default rates, recovery rates or exposures, fluctuate 
substantially in time even on a high portfolio aggregation level jj- This results in 
large unexpected losses in loan portfolios, for which banks are required to hold equity 
capital as a loss buffer. To determine the appropriate level of equity capital for banks' 
loan portfolios is one main focus of the regulatory consultive process known as Basel 
II [2]. Accordingly, credit risk modelling has been a focus of intense research in recent 
years [3 II © EJ Ell El H although considering 

the risk premium when pricing interest rates goes back some time |22j . 

The assessment of credit ratings by assessment agencies such as Moody's and 
S&P allow some statistical assessment of the credit quality of individual offerings or 
particular companies. However, it is clearly essential when considering the risk of a 
basket of loans that the correlations between the members of the portfolio are taken 
into account. One systematic approach is to replace the number of firms in a portfolio 
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with an effective number of independent firms [21]. By boosting the contribution of 
each firm to keep the mean loss constant, this introduces a larger variance of losses, 
in at attempt to capture the risk caused by correlations between firms. JPMorgan's 
Credit Metrics approach 3j (see also Credit Suisse Financial Products's CreditRisk 4 " [I] 
and jinj for a detailed comparison of the two) tries to model the correlations between 
firms in credit quality using the observable correlations in equity value of the firms. 
An intuitively appealing approach is to assume that the default intensity depends on 
some set of macroscopic economic factors (e.g. interest rates, growth rates, oil prices 
etc.), the so-called reduced form model [7||Hj. Thus the default rates of different firms 
are coupled via some limited number of factors, but given the factors the default rates 
are independent. In structural models [221 123] - the dependence on macro-economic 
factors is understood in terms of correlations in the dynamics of asset returns of 
different companies, leading to correlations in default rates via correlated dynamics 
of returns. More involved approaches have modelled interactions between firms in the 
wider economy by introducing changes to a firm's default intensity upon the default of 
another firm [S], or via a copula function [10j . The quantification of these correlations 
via simulations is discussed in |14j . 

The main purpose of the present contribution is specifically to expand upon 
recent modelling and analytic descriptions of the influence of counter-party risk 
[HlliamElQSllIHlimilBlEI]. Counter-party risk addresses the fact that a given 
firm's economic health is strongly influenced by the performance within the network of 
companies with which it has direct economic interactions. Economic interactions are 
here to be read in a broad sense as any form of relation which is likely to mutually 
influence performance, though not necessarily in symmetric ways. Specifically, a 
defaulting firm within this network of counter-parties will affect a company's own default 
probability — reducing it, if the defaulting firm was a competitor, or increasing it, if 
the relation was of a cooperative nature. When the default probability is increased, this 
process is known as credit contagion and has been considered in e.g. fHJ ^] while the 
incorporation of counter-party risk into a reduced form model was introduced in [0]. In 
what follows we consider the dynamics of individual firm defaults and their influence on 
loss distributions. More subtle effects such as credit quality migration are, as yet, not 
taken into account. 

The importance of direct functional interactions in the analysis of risk is not 
restricted to credit risk. In fact the role of interactions is much more obvious in the 
context of operational risk, where sequential, functionally induced failures of mutually 
dependent processes constitute one of the main sources for operational risk. Indeed, an 
attempt to explore the consequences of interactions for quantifying the capital buffer 
necessary to cover operational risk j2Hj has provided major ingredients for the approach 
to credit risk modelling started by Neu and one of us |2*T] . 

In the present paper we provide an analytic solution to the dynamical description 
of counter-party risk within a heterogeneous, functionally defined network of interacting 
firms (to be referred to as economy in what follows), see e.g. [TH], in the spirit of |21j . 
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We generalize the analysis of that study to capture effects of cooperative as well as 
competitive business relations within the economy, and we solve the model for a wide 
degree of dilution of the network of economic dependencies in the sense that we assume 
each company in the net to have business relations only with a (randomly chosen) subset 
of the full set of companies. In the present investigation we will always consider the case 
where the number of interaction partners of each company is large, and for simplicity we 
shall restrict ourselves here to the case where the graph defining economic connectivity is 
a Poisson degree distributed Erdos-Renyi random graph [2E]. More realistic connectivity 
distributions reflecting the different connectivity patterns of large and small players in 
an economy, taking into account small-world effects and fat tailed degree distributions 
[27] can be handled by methods similar to those used in the present investigation [2Hll29| . 
but will be studied in a separate paper. 

As in [21], the model parameters are unconditional and conditional default 
probabilities, which may be thought of as being obtained via a suitable rating procedure. 
One of the virtues of the present analytic investigation is to highlight the fact that the 
collective behaviour of the system, which ultimately determines the the loss-distribution 
on an economy-wide scale, is fairly insensitive to detail. That is, it does not depend on 
getting individual dependencies correct, but only on the overall distribution of of the 
unconditional and conditional default probabilities. Our main result is to demonstrate 
that the effects of economic interactions — while relatively weak in typical scenarios - 
are pronounced in situations of economic stress, and thus lead to a substantial fattening 
of the tails of loss distributions even in large loan portfolios. 

The remainder of the paper is organized as follows. In Sec. El we define our model, 
and specify the stochastic setting for our analytic investigation. The relation between 
the model parameters and conditional and unconditional default probabilities as used 
in [21] is briefly reviewed to make the paper self-contained. Sec. H3 describes a heuristic 
solution for the dynamical evolution of the fraction of defaulted companies over a risk 
horizon of one year, starting from an overall healthy situation, a scenario appropriate for 
the analysis of credit risk. A formal solution in terms of a generating function approach 
(GFA) [30 , which provides a full justification for the heuristic solution is relegated 
to |Appendix A Both solution methods are based upon techniques developed in the 
statistical mechanical analysis of dilute neural networks [21]. In Sec. E]a phase diagram 
distinguishing regions in parameter space in which economic interactions can lead to a 
collective acceleration of the economy-wide default rate in situations of economic stress 
from regions where such acceleration is impossible is computed. Distributions of annual 
fractions of defaulted companies as well as loss distributions, both economy-wide and for 
finite loan-portfolios are also computed and compared with simulations. Sec. 03 finally 
summarizes our findings and discusses their main implications for the analysis of credit 
risk 
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2. Model Definitions 

In this section we define a statistical model that attempts to capture the effects 
of counter-party risk on credit contagion. In contrast to approaches based on 
microeconomics, and in keeping with the framework discussed in the introduction, we 
allow the firms' wealth, macro-economic factors and interactions between firms to all 
be described probabilistically. This is due to our focus on the characteristic change in 
behaviour caused by examining interactions between firms in the wider economy. 

We analyse an economy which consists of N firms. The state of each firm j at a 
given time t is described by its 'wealth' Wu which we take to be the difference between 
its assets and its liabilities. Accordingly, a company defaults, if its wealth Wu falls 
below zero. As indicated in the introduction, we are interested in the influence of 
economic partners on firms' own performance, specifically on their default probability. 
We understand the notions of firm and economic partner in a very wide sense: a firm 
could be any economic entity, a manufacturer, a service provider, a trader of goods or 
services, etc. It could also be an individual (we shall often use the generic term 'node' 
to designate these entities). Two firms are partners if the state of one has a material 
effect upon the other, e.g. one is the supplier of the other, performs outsourced services, 
there exists substantial loans or other financial commitments between the two, or they 
compete in the same market. For instance, if a major manufacturer of PCs were to go 
out of business tomorrow, this would inevitably have a material impact on the economic 
performance (i.e. wealth) of other companies operating in the same industry sector. 
It would on the one hand lead to a deterioration in the financial viability of most of 
the suppliers or service providers of the PC manufacturer in question — including in 
particular its own work force! — but on the other hand it would improve the situation 
for competing producers of PCs, in that they could profit by taking over a share of the 
defaulted company's market. 

We are interested in quantifying the effect of these interactions from the perspective 
of a lending bank which would be required to set aside a sufficient amount of capital 
to cover losses incurred by defaults of its obligors. Another perspective might be that 
of a central bank, which would base monetary policy decisions in part on their impact 
on expected default rates at an economy wide scale. The typical risk horizon in these 
contexts would be one year. 

For simplicity, we assume that that within the risk horizon of one year node 
% experiences an interaction-induced material change of its wealth only if one of its 
business partners, say j, defaults. In order to formalize this in a dynamical description, 
we introduce binary indicator variable rijt which indicates whether node j is solvent at 
time t {rijt = 0) or has defaulted (rijt = 1). 

The value of the ith node's wealth at time t, W it , is thus taken to be of the form 

n 

W it = JijUjt -r}u (1) 

Here J^- denotes the change in i's wealth which would be induced by a default of node 
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j. One would have > if j is a cooperative partner of i, whereas J^- < if j is a 
competitor, while Jjj = if there is no direct influence of j on z. By $j we designate z's 
initial wealth at the beginning of the year, and rju are random fluctuations caused by 
both, external macro-economic factors (an expanding/shrinking economy, an oil price 
spike, market sentiment etc), and firm specific actions or events. 

We take the initial state of the economy to be a set of solvent firms, n iQ = for all 
i, (one could view this as a definition) and say that a firm defaults at time t if Wu < 0. 
We define our dynamics such that if a firm goes bankrupt it does not recover within a 
risk horizon of one year, so the bankrupt state is absorbing. Thus the dynamics of the 
firms state is given by the equation 

n it +i = n it + (1 - n it )Q ^ ^ J i:j n jt - + 77^ (2) 

where ©(. . .) is the Heavyside function. The time step in this dynamical rule will be 
taken to represent one month. 

We choose the r\a to be Gaussian distributed. Without loss of generality they can 
- by suitably rescaling the and the — be chosen to have unit variance. We follow 
widespread practice IHj to account for common fluctuating macro-economic factors 
by choosing the rju to be correlated for different i. This could be achieved by taking r/u 
to be of the form 

K 

Vit = O-iiu + PikYkt (3) 
k=l 

with uncorrelated Gaussian unit variance white noises and {1^}, the former 
describing firm specific wealth fluctuations, whereas the latter could account for the 
relative effects of fluctuations common to industry sectors, regions, or countries, with 
prefactors Oi and (3^ describing the relative importance of these fluctuations on i. In 
what follows we restrict ourselves to a minimal variant of this set-up by finally choosing 



Vit = y/p Vot + V 1 -P tit ■ ( 4 ) 
We shall simplify matters further by assuming that the common economic factor r] ot is 
slow and take it to be constant 7]o t = 770 within a risk horizon of one year. One-factor 
models of this type also feature in the regulatory framework laid out in the Basel II 
accord |2j. 

None of the simplifying assumptions are necessary for our analysis to go through; 
the generating function formalism given in |Appendix A in particular can easily handle 



more general cases. However, the simplified setting is sufficient to highlight the 
important effects of interactions on credit risk, and it does lead to a greatly simplified 
macroscopic description of the system, as we will see in Sec. El 

As for the which describe the loss or gain of node i due to a default of node j, 
in the present paper we will investigate them in a probabilistic setting. It will be useful 
to disentangle the presence or absence of an interaction from its strength by writing 

Jij — c ijJij (5) 
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where c y - G {0, 1} describes the absence or presence of a connection j — > i, while 
describes its magnitude, both of which are assumed to be fixed. It is reasonable to 
assume that connectivity is a symmetric relation, cy- = Cji, whereas there is no reason 
to suppose symmetry of the magnitudes of mutual influences. Considering the case of 
a small supplier with one large company taking the majority of its orders, if the larger 
company defaults then the small supplier may well go bust too. However, if the small 
supplier defaults then the large company is less likely to suffer terminal financial distress, 
so in general jy ^ Jji. 

Specifically, we assume a random connectivity pattern described by 

p ( c ij) = J^ 6 ^,i + (l - <Jcy,o , i < j , = c^ (6) 

and we will be interested in the limit of a large economy (the thermodynamic limit), in 
which the average connectivity c of each node is itself large, N — > oo, c — > oo. We will 
initially be concerned with the extremely diluted regime, where c/N — > 0, taking e.g. 
c = 0(log(iV)). These assumptions have important consequences for the structure of 
the graph defining the connectivities, namely that each node feels the effects of a large 
number of other nodes (so that limit theorems will allow us to describe the overall effects 
of interactions) and, for the extremely diluted regime, that there are only a finite number 
of loops of finite length even in the infinite economy limit. The graph of interactions 
between companies for finite N is just an Erdos-Renyi random graph [26 . 

The magnitudes Jy of the interactions will also be taken as fixed random quantities. 
In order to allow the taking of the thermodynamic limit as described, the mean and 
fluctuations of the Jy must scale in a suitable way with the connectivity c. Quite 
generally, we must have 

Jij = — + ~^ X ij (?) 

in which the Xij are zero-mean unit-variance random variables. The scaling of mean 
and variance of the Jy is given by the parameters Jo and J respectively. If Jo > there 
will be a net cooperative tendency within the economy, which seems to be a reasonable 
assumption. Finally, we will need to assume that all moments of the are finite and 
we will choose the to be independent in pairs 

Xij = , xfj = 1 , XijXji = a , XijXki = otherwise . (8) 

The parameter a, (— 1 < a < 1) describes the degree of correlations between J y - and 
Jji. Strictly symmetric interactions are obtained only for a = 1 

At this point let us briefly recall that, after rescaling as described, the model 
parameters i?j and Jy have a clear meaning in terms of unconditional and conditional 
default probabilities [21]. We denote by n t the values of all indicator variables in the 
economy at time t, assuming riu = 0. Then by integrating over the unit variance 
Gaussian rju in (J2J one obtains the conditional probability for node i to default within 
a month given a configuration n t of non-defaulted and defaulted firms in the economy 
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at time t as, 



Prob (n 



it+1 



n = ^{J2 j ij n jt ~ ^ 



with = |[1 + erf(x/-\/2)] the cumulative normal distribution. Thus the 

unconditional probability pi of default of i within a month in an otherwise healthy 
economy and the conditional probability pi\j for a default of i within a month, given j 
and only j has defaulted before are given by 

Pi = Prob (^t+i = l| fat = 0}) = $(-#i) , (9) 

Piy = Prob = 1 rijt = 1, {rife(^-)t = 0}) = Q^Jij - . (10) 

These relations may be inverted to express the model parameters in terms of conditional 
and unconditional default probabilities — quantities that would be estimated in a rating 
procedure — as 

■&i = -<S>-\pi) , Ja = Q- 1 ^) - ^{Pi) ■ (11) 

While characterising the default of companies is of interest, our primary concern 
is to examine the distribution of losses accrued over our one year time frame, both in 
the economy at large, and in a portfolio made up of a finite number of firms within the 
economy. We assume that the losses caused by default are independent of the month of 
default, and then examine two different cases. The first simpler case is that the losses at 
firm i, given that firm i defaults, are uncorrelated with any other variables. The second, 
perhaps more interesting case, is that the losses at firm i are random but are correlated 
with the initial monetary reserves The intuitive reasons are that if a firm has more 
cash, then the default is less anticipated, and thus will be less priced in by the market; 
or the firm has larger credit lines and so will default on a larger amount; and finally the 
firm is likely to be larger, and hence cause a larger loss. 

3. Heuristic Solution 

In the present section we show that our model has a relatively simple solution that can be 
obtained by qualitative probabilistic reasoning, appealing to statistical limit theorems. 
This solution turns out to be exact, as we show using a more involved generating function 
formalism in Appendix A Both types of argument have been developed in the analysis 



of the statistical mechanics of disordered systems, and in particular neural network 
models while for a more general introduction to emergent collective behaviour see 
e.g. [^21. R eca H the microscopic dynamics as defined by ©. The complications are 
due to the interactions between firms, namely that the state of a given firm i at time t, 
depends on the state of the neighbours of i for times t' < t which in turn depend on i at 
times t" <t'<t. In general this feedback prohibits straightforward analysis, and indeed, 
it led Jarrow and Yu jH] to eliminate this feedback explicitly by considering an economy 
of two types of firms: primary firms whose default depended only on macro-economic 
factors and secondary firms whose default depended on macro-economic factors and the 
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default of primary firms. However, due to the specific structure of our model we are 
able to push the analysis further. By definition, the overall effect of interaction terms 
on company i at time t is given by the local field hu = J2j Jijn-jt- From the statistics of 
the interactions J^- given by ©-© we see that each firm i is connected to, on average, 
c other firms. Since we consider the large c limit, this means that we could evaluate the 
statistics of h it by appeal to the law of large numbers and the central limit theorem if 
the contributions to h it were independent, or at least sufficiently weakly correlated. 

At first sight we cannot expect this condition to hold if we have some degree of 
symmetry in the interactions, i.e. for a^O, even in the extremely diluted regime. Note 
that there are two ways in which the njt of the neighbours interacting with i may become 
correlated through the dynamics: either they influence each other through firm i, or not 
through firm % but though some loop of interactions in the economy. In the extremely 
diluted regime, correlations between the neighbours j of % cannot build up in finite time 
(within the risk horizon) via loops not involving i, since due to the scaling almost all 
loops are very long. With symmetry in the interactions, correlations between the nj t 
could in principle be induced by the dynamics of i. However, as long as = 0, the 
njt clearly cannot influence each other through site i, whereas once nu = 1, then firm % 
is in the absorbing state, and correlations it induces on the dynamics of its neighbours, 
have become irrelevant for its own microscopic dynamics (0). Thus limit theorems can 
be used after all to solve the macroscopic dynamics of the system, despite a possible 
symmetry in the interactions. 

Returning to the dynamical evolution equation (j2J, we observe that the coupling 
of a node to the economy is via the local field 

Kt = X! Jij n 3t — ~~ Y C ij n jt + —j= Y C ij X ij n jt ) (12) 
3 ° 3 VC j 

which is a sum of random quantities (with randomness both due to the Gaussian 
fluctuating forces (the {r] it }, respectively the {£it}), and due to the heterogeneity of 
the environment). The first contribution is a sum of terms of non- vanishing average. 
By the law of large numbers this sum converges to the sum of averages in the large c 
limit, 

Jo Jo v-» j r Jo v-» ~, r T 1 



h it = — £ °i3 n 3t ^—Y, C i3i n 3t) ~ — £ °ij ( U 3t) = J0J7 J2 M 
3 3 3 3 

in which angled brackets (...) denote an average over the fluctuating forces, and the 
overbar (. . .) an average over the Jy, i.e., the and the Xij. An approximation is 
made by assuming negligible correlations between the Cjj and the (nj t ) induced by the 
heterogeneity of the interactions. The second contribution to (|12|) is a sum of random 
variables with zero mean, which we have argued are sufficiently weakly correlated for 
the central limit theorem to apply for describing the statistics of their sum. Thus the 
sum ^ 
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J' 



J 2 



{{Shu) 2 ) — /] CijCik x ii x ik( n it' r lkt) — 22 CijCik%ij%ik( n jt n kt) 



An approximation based on assuming negligible correlations has been made as for the 
first contributions. Thus the local field h it is a Gaussian with mean h® t and variance 
((5h it ) 2 ) both scaling with the average fraction of defaulted nodes in the economy. By the 
law of large numbers this average fraction will be typically realized in a large economy, 
i.e. we have 



m t 



1 

N 



E 



n 



.it 



N 



J2 M 



(13) 



3 

in the large N limit. The dynamics of the fraction of defaulted nodes then follows from 



m *+i = y E = m * + Jj E ( X - nit)® (ha ~ &i + y/PVo + \A ~ P & 



it j 



(14) 



where the one factor noise model (J1J) has been used. 

The sum in (JT3j) is evaluated as a sum of averages over joint n it , ha, and 
distribution by the law of large numbers. We exploit the fact that n^, and /i^ 
are uncorrelated. Noting that the sum h it + a/1 — p is Gaussian with mean Jo^i and 
variance 1 — p + J 2 m t , and taking into account that n^-averages, depend on i through 
(nit) = (nt)(#i), we find 



1 _1 

m t+ i = m t + — 2^ - 



(n t ) 



1 + erf 



^2(1 -p+J 2 m t ) 



This version can be understood as an average over the d distribution 

1 



N 



J2W-$i) 



which maps onto a distribution of unconditional default probabilities as discussed above. 
Denoting that average by (. . we finally get the following evolution equation for the 
macroscopic fraction of defaulted companies in the economy 



m t +i =m t + 



1 " (nt) 



erf 



Jam + y/p Vo - $ \ 



(15) 



'2(1 -p + J 2 m t ) J 

We have thus an explicit dynamic equation for the macroscopic fraction of defaulted 
nodes in the economy. It involves first propagating ^-dependent default probabilities 



via 



( n t)w 



1 + erf 




J m t + y/p f] - g 
- p + J 2 m t ) _ 



(16) 
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which depends only on m t} thereafter performing an integral over the $ distribution to 
obtain the updated fraction m t+ i of defaulted nodes given in (fI5jl. 

The heuristic solution of the macroscopic dynamics 1)1 5|) . (jlBj) presented here is 
based on independence assumptions which are not easily justified in a rigorous way via 
the probabilistic reasoning presented above. However, the solution is supported in full 
detail by an exact analysis based on generating functions presented in |Appendix A 



4. Results 



In the present section we explore the consequences of our theory. We studied the 
dynamics and computed loss distributions for an economy in which the parameters $j 
determining unconditional monthly default probabilities according to (jHJ) are normally 
distributed with mean $o — 3, and variance a\ = 0.01 so that typical monthly 
default probabilities are in the 5 xlO -4 range. Except when stated otherwise we shall 
use p = 0.15 for the parameter describing the relative importance of economy- wide 
fluctuations, a value that is considered to be in an economically acceptable range. 

In Fig la we show the evolution of the typical fraction of defaulted firms over a 
risk horizon of 12 months for various settings of the interaction parameters Jo and J; 
the typical fraction is computed by choosing the most-probable value r] = for the 
economy-wide influence on the dynamics. Fig lb shows the probability density of the 
end of year fraction of defaulted firms driven by fluctuations in economic conditions. It 
is obvious that interactions cause a significant fattening of the tail of the density at large 
values of this fraction, which is a clear indication of the significance of counter-party 
risk in particular in situations of economic stress. 




Figure 1. Left: Typical fraction of defaulted companies as a function of time for 
(J , J) = (0,0), (1,0), (0, 1), and (1,1) (bottom to top). Right: Distribution of the 
fraction of defaulted companies at t = 12 for (Jo, J) = (0,0), and (1,1). Analytic 
curves are overlaid with results of a simulation in which the distribution is obtained 
by computing the fraction of defaulted companies for randomly sampled 770. 



In order to assess whether interactions can lead to a collective acceleration of the 
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rate of defaults we look at the discrete second derivatives 

A t = m t+1 + m t _i - 2m t 

which are always negative for the non-interacting system and maximal at t = 1, 
irrespectively of 770- Interactions can lead to a collective acceleration of the rate of 
defaults signified by the possibility that the A t may become positive in unfavourable 
economic conditions. We define the region in parameter space in which collective 
acceleration of default rates can occur by the condition that 

Prob{A t > 0} > (17) 

for some t, with 1 < t < 11. The concavity of the error function for positive arguments 
entails that effects of collective acceleration are always strongest at t — 1. Evaluating 
this condition for various values of the parameter p describing the coupling to the overall 
economy, we get lines shown in the phase diagram Fig 2a. Note that the influence of p 
is very weak in interesting region of low p values. Note also that Ai values are typically 
positive but very small in the region near the phase boundaries shown, as illustrated in 
In Fig 2b where we exhibit the distribution of discrete second derivatives at t — 1, just 
inside the phase where acceleration of default rates is observed. The tail of negative A x 
is found to extend to significantly larger values. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 -0.0005 -0.0004 -0.0003 -0.0002 -0.0001 0.0001 0.0002 



J 4, 

Figure 2. Left: Phase boundaries separating regions without collective acceleration of 
default rates from regions where acceleration occurs, for p = 0.15, 0.3 and 0.8 (bottom 
to top). Distribution of discrete second derivative Ai = to 2 + too — 2m i just within 
the phase with accelerating default rates for ( J , J) = (0.5, 0.5) (right). 

The quantity of central importance from the point of view of credit risk analysis is 
of course the distribution of losses. Let ^ denote the loss that would be incurred by a 
default of node i. Then the loss per node for a given state r/ of the economy is 

HVo) = ^J2 n ^ (18) 

i 

where £j is randomly sampled from the loss distribution for node i. We assume that 
the £i are independent of the stochastic evolution. In the large system limit, the loss 
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per node at given value of r] describing the influence of the overall economy is a non- 
fluctuating quantity, as it is itself an (empirical) average taken over an (infinitely) large 

system, 

L( Vo ) = lim ^EMW^i) = / dM^Mw'm (19) 

N—>oo iv • J 
i 

by the law of large numbers, where £ = is the mean of the loss distribution for 
a node. If loss distributions were identical for each node, with means independent of 
default probabilities = £o, then the distribution of losses driven by the fluctuations 
of the economic stresses would simply replicate the distribution of the fraction of 
defaulted firms. 

The situation is different if loss-distributions are correlated with default 
probabilities. As an example we consider the case where average losses are inversely 
proportional to the unconditional default probabilities Pd($i) = Pi introduced in (JHJ). 

£(&) = ^— (20) 

with a parameter e introduced as a regularizer to prevent divergence as pd — > 
0. That is, the contribution to the total losses incurred by defaulting firms with 
different unconditional default probabilities is approximately uniform over the default 
probabilities. In our model we have 

p„(0) = i[l-erf(0/V2)] (21) 

Fig. 3 shows loss distributions for such a situation. The analytic curves are computed 
by noting that the losses per node are monotone increasing functions of rjo which is itself 
jV(0, 1). Integrated loss distributions are thus simply obtained using error functions 

Prob[LM <L] = i[l + erf(r /0 (L)/v / 2)] 

where rj (L) is the 7/o-value giving rise to loss L per node. The probability density 
function is obtained via a single numerical differentiation. 

It would be of some interest to know whether finite sample fluctuations could 
possibly upset the picture seen so far. To study this issue we look at the losses per 
node of a finite sample randomly drawn from the nodes of a large economy, 

i M 

M») = i7E n u(Wi) ( 22 ) 



Writing this as 



M i=1 



I M 



M»3b) = TpE(< n i3>wW 

1V1 i=l 



i=i 

i M 



L(Vo) + L 2 ( Vo , {■&<}) + L 3 ( Vo , {^}, Wit}) (23) 
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Figure 3. Loss-distribution per node for the infinite system with = l/(e +Pd($)) 
at £ = 0.005. Analytic curves are overlaid with simulation results. Lower curve: 
non-interacting system, upper curve: interacting system with (Jo, J) = (1, 1). 



we see that it has three components. The first is simply the expectation value describing 
the loss per node at given t]q in an infinite system, the third, L3 has zero mean and is 
expected to be Gaussian at large M, with variance scaling as M -1 describing the noise 
induced fluctuations about the average for a given collection of {$i}, while the second, 
L2 — also a a zero mean Gaussian of variance scaling as M _1 at large M — describes 
the finite sample fluctuations of this average. Since the collection of {^j is fixed, 
these Gaussians are correlated for different r]o . While an analytic evaluation of the loss 
distribution may still be feasible in principle, it would become very involved in practice. 

An approximation to the finite size computation is obtained by assuming that losses 
per node at given rj are Normally distributed about their infinite system 770-dependent 
mean with [combining L 2 and L 3 ] variance (also 7/0-dependent) 

4/ = M- 1 (({n 12 ) m ¥W)), - (Mi»)>^)^ , (24) 

which is an annealed approximation which ignores that the parameters of the individual 
loss distributions of the nodes in question remain fixed. The results of an evaluation 
along this line are shown in Fig 4, using the scaling ()2U|) of average losses used above; 
the approximation suggests that finite size fluctuations give rise to fatter tails in the 
loss distributions, though the effect is negligible in the interacting system except at the 
extreme end of the loss distribution, and small but a bit more pronounced in the non- 
interacting system. A comparison with a simulation shows that the effects of fluctuations 
are slightly underestimated in our approximation. Note, however, that loan portfolios 
of typical banks usually contain orders of magnitude more debtors than the M = 100 
considered in the present example. 

Let us finally look at the so-called Value at Risk in terms of which the capital buffer 
that banks are required to hold to cover risk is often expressed. It is defined as 

VaR g = (Q q [L] - (L)) e~ rT (25) 
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Figure 4. Loss-distribution per node in the large system system limit with = 
l/(e + Pd($)) at e = 0.005 and for a finite sample of M = 100 companies randomly 
taken from the ensemble. For the finite sample, individual loss distributions are taken 
to be flat in the range [0, 2l{$\. Right: The same analytic curves, overlaid with a 
simulation result for the finite sample. Lower curves correspond in both cases to the 
non-interacting system, upper curves to the interacting system with (Jo, J) — (1, 1). 

in which Q g [X] is the g-quantile of the loss distribution at time T, i.e. the loss that is 
not exceeded with probability q 

Prob(L < Q q [L\) = q , 

while (L) is the average loss, and r denotes a risk free interest rate. To highlight 
the effects introduced by economic interactions we take the ratio of the value at risk 
computed for economies with and without functional economic interactions VaR/VaRo, 
both computed at confidence level q = 0.999 as required by the Basel II regulations 
Taking this ratio also eliminates the dependence on the interest rate r. 

In Fig 5 we display this ratio alongside with analogous ratios of average losses 
(L)/(L) for comparison. 

As perhaps may be anticipated in view of results displayed in Figs 3 and 4, the 
Value at Risk is significantly more sensitive to functional interactions in an economy, 
than the average losses are. This is understandable, as VaR probes the tails of loss 
distributions, while average losses will be determined mostly by typical results. 

For the results displayed in Fig 5, unconditional default probabilities where not 
adjusted with the strength of the interactions so as to keep the average annual default 
probability constant. In Fig 6, therefore, we take this extra step, displaying an 
analogous ratio of the Value at Risk of interacting and non-interacting economies, 
where now unconditional default probabilities in the interacting system are adjusted 
in such a way that the average (interaction-renormalized) default probability stays 
constant — at the level chosen for the non-interacting system. To keep matters simple, 
a homogeneous portfolio, with firm-independent unconditional default probabilities 
and firm independent average losses was chosen. Clearly the interaction-induced 
enhancement of the Value at Risk is rather close to the corresponding enhancement 
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R 

Figure 5. Ratio of Value at Risk for systems with and without functional interaction 
as a function of the strength of the interaction (upper curves) . The analogous ratio 
for average losses is also shown in each case (lower curves). The curves shown in the 
three figures are evaluated along straight lines in the Jo — J plane, and the parameter 
R measures a distance from the origin R = y/ Jq + J 2 , The three figure correspond to 
the lines Jo — (upper left), J — (upper right) and J /J = 1 (lower). 

computed without adjustment of the unconditional default probabilities. 

To summarize, the capital buffer that banks are required to hold according to 
the Basel II regulations [2j to cover credit risk is significantly underestimated when 
interaction effects in an economy are not taken into account. It is important to note 
that this is true already in the regime in which interactions are too weak to cause an 
overall acceleration of default rates, as can be seen by comparing the phase diagram 
Fig. 2a with results for the Value at Risk displayed in Figs. 5 and 6. 

5. Conclusion 

In conclusion, we have studied the effects of economic interactions on credit risks. 
Though non-equilibrium initial conditions and the fact that the credit-risk problem 
has an absorbing state would at first sight appear to complicate the analysis, we found, 
quite to our own surprise, that in particular the presence of the absorbing state simplifies 
the analysis considerably, as it removes the non-Markovian effects in the macroscopic 
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Figure 6. Value at Risk relative to Value at Risk in systems without functional 
interaction as a function of the strength R = y/ Jq + J 2 of the interaction for Jo/ J = 1 , 
with unconditional default probabilities in the interacting system adjusted as a function 
of the strength R of the interaction to keep the annual average default probability 
constant (in the present case at a value close to 3%.). As in Fig. 5, the VaR ratio is 
computed at q — 0.999. 

dynamics that would otherwise be present in systems with some degree of symmetry in 
the interactions. While the limit of extreme dilution simplified the reasoning within the 
heuristic solution, we saw in the generating function analysis that the assumption of 
extreme dilution could be dispensed with. So although the rather heavy machinery of 
non-equilibrium disordered systems theory is required to rigorously treat the model (due 
to asymmetry in the inter- firm dependencies and the initial conditions), the resulting 
effective single-firm process is remarkably simple. This has obvious practical benefits in 
terms of computational efficiency 

We have seen that the effects of economic interactions are relatively weak in typical 
economic scenarios, but they are pronounced in situations of economic stress, and thus 
lead to a substantial fattening of the tails of loss distributions in large loan portfolios. 
This leads to significant increases in the Value at Risk, i.e. the capital that must be 
held as a loss buffer, when compared to the non-interacting theory. Importantly, this 
conclusion remains valid even in the case where there is no overall acceleration in default 
rates, c.f. Fig. 2a and Figs. 5 and 6. 

It is worth paraphrasing these last observations as they address a point of key 
importance. While credit risk models that do not take direct economic interactions 
into account can provide a very reasonable fit, when calibrated on historical data which 
reflect normal economic conditions, their predictions would be entirely inadequate when 
it comes to estimating default rates and losses in situations of significant economic stress. 

Note that the model presented here is suitable for detailed and comprehensive stress 
testing, as explicitly demanded within the regulatory framework of the Basel II accord 

The issue of stress testing was addressed in greater detail when the present model 
was first introduced in [2*T] . 

The patterns of economic interactions studied in the present paper are described 
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by an Erdos-Renyi random graph. The large connectivity limit considered in the 
present investigation further entails that there is no pronounced heterogeneity in the 
sets of economic partners of any one given node. Connectivity distributions other than 
Poisson can, however, be handled by suitably adapting the generating function approach 
explained in |Appendix A along the lines developed in [2H] , and will be investigated in 



a separate publication 33J. In terms of model fitting there appear to be a vast number 
of free parameters in terms of the interactions {Jij} between firms. However, it is 
important to realise that to understand the macroscopic behaviour, here only their low 
order statistics are relevant. 

In the present investigation, we restricted ourselves to analysing the effects of 
interactions on default-dynamics and, via default rates, on loss distributions. More 
subtle effects such as credit-quality migration are, as yet, not taken into account, but 
could be modelled along similar lines using the dynamics of interacting multi-state 
indicator variables. Further assumptions concerning details of such models would be 
required, however, and the full complexity of non-Markovian dynamics would resurface 
in such an analysis. 

Appendix A. Generating Function Analysis 

In this appendix we describe the generating function approach (GFA) to solve our model, 
giving full justification to the arguments used in section El The reasoning is relatively 
standard; we include it here to make the paper reasonably self-contained 

Appendix A.l. The Generating Function for Correlation Functions 

First we introduce the generating function at fixed value of the macro-economic force 

Z[^\r] } = ^SwEi^^ (A.l) 
where the angled brackets denote averages over the microscopic dynamics (J2J) of rii, i.e. 

^[VM= E i 5 [n ,...,n 12 ]e- i S"oS*^ n » , (A.2) 

n ,...,ni2 

with P[no, . . • , TI12] denoting the probability of a sequence of configurations of the entire 
set of interacting firms over the risk period of 12 months. The generating function can 
be used to compute expectation values and correlation functions via differentiations 
with respect to the source fields ipu, 



(jHt) = 1 



dipi, 



, {riisnjt} = 1 — 

^0 d^ is d^ jt 



It is expected that correlation functions averaged over the randomness in the couplings 
Jij are dominated by typical realizations of the disorder, hence to describe typical results 
an average of the generating function over the disorder, 

ZWM = J II dP (% J ji) Z[1>\Vo) (A.3) 

i<j 
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is computed. 

To proceed, the path-probability P[n Q , . . . , n 12 ] at given t)q is expressed in terms of 
transition probabilities of the Markovian dynamics, 



11 



P[n Q , n X2 ] = P(n ) ]J P(n t +i\n t ) 



t=o 



where 



with 



ft = n it + (1 - n it )Q J2 J ij n Jt + V 1 _ P &t + VP Vo ~ & 



(A.4) 



(A.5) 



The £it integrations appearing in the transition probabilities and the average over the 
Jij distribution are facilitated by utilizing ^-distributions to 'extract' the £ it and the Jy 
from the Heaviside function in f it , using 



1 = J du it 5 \u it - Jij n jt - Y 1 - P C 

This gives 

P(n t+ i\n t/ 



e V 3 

2vr 



jt- 



J -q du it du it ~ l i Eil2 it~' lil it( u it-Y.j J io n jt)\ -q^ ^ (A 6) 



with now 

fit = n it + (1 - n it )e (u it + sfp rj - tf;) 
Inserting into the generating function, we get 



(A.7) 



Z[#fc]= £ ^K)/n^^exp^ 

no, ...,ni2 J it Z7F 



\ m it) 



H 3n it+1 ,f it 



The disorder average affects the J^; it factorises in the pairs and involves the term 



II Aj = II ex P V H (u it Ji^t + UjtJjiUitj 

(i,j) i<j I * 



Here the superscripts c and x indicate averages over the Cjj and the in the 
according to the statistics Performing the Qj average, one gets 



n Aj = n j 1 + ^ ex p { + s 



£ru.,tn if j - 1 
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The exponential is expanded using c 3> 1. Using (jSJ), keeping dominant terms and 
re-exponentiating the result one obtains 



n 



exp < N 



J Y hm t + — Y [QstQst + aG st G ts 

t ^ s,t 



which depends only on the macro-variables 



N 



Y™u , rn t 



1 

N 



Y n it 



Qst = Y ' l UisiUit , q st = — Y n is n it , Gst 



Y, ^is n it 



(A.8) 
(A.9) 



We thus have 



Z[*o] 



E ^("o)/n^expfe 



no,...,n 12 



— — [lU it ) - lUuUit - liPitTlit 



+N 



J Y k t m t + ~W Y [QstQst + aG st Gts 



n<Wij 



Site factorisation in Z^lr/o] is achieved as usual by writing it as an integral over the 
macro-variables, using 5-function identities of the form 



1 



d(Nm t ) 5^Nm t - Y 



n 



.it 



dm t drht im 
e 



2it/N 



and analogous ones for the k t , q st , Q st , and the G st to compute densities of state. This 
results in the following compact expression for the average generating function 



= Jv{.. .} exp {iV[$ + * + S]} 



(A.10) 



in which T>{. . .} stands for an differentials of all order parameters introduced in flA.8|) . 
flA.9|) and their conjugate (hatted) parameters introduced via Fourier-representations of 
(^-functions. The functions $, and S appearing in (jA.10|) are given by 



$ = Jo Y k t m t + ~W Y [QstQst + aG st G ts 

t ^ s,t 

* = lYi^tmt + k t k t ] + iY^QstQst + QstQst + G st G st ] 



(A.ll) 
(A.12) 



St 



t?e io «e /n 



N 



dutdut 



2tt 



exp - S - iY^itrit) n<Wi,/« -(A. 13) 



with S denoting the 'dynamic action' 



5 = E 



1 - P t;.~. \2 



m t ) +iu t u t +imtnt+iktiut 



St 



q st n s n t +Q s tiu s iut + G st iu s nt 



(A. 14) 



The third contribution, H, in (|A.10|) describes an ensemble of independent single site 
dynamical problems. Thus, to leading order in N we have written our generating 
function in terms of an integral which may be computed via a saddle point argument. 
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im t = 


- Joh 


ik t = 


-J m t 




J 2 


iQst 


J 2 


iq s t = 


T^Qst 


— — 2~^ s * s * 


m t = 


i 


k t = 


i 


q s t = 


i 


Qst - 


i 


G s t — 


•i 







Appendix A. 2. Saddle Point Problem 

In the saddle-point, variation of our observables gives 

-aJ 2 G ts (A. 15) 

(A.16) 
(A.17) 
(A.18) 

with (•••)(«) denoting averages evaluated wrt effective single site dynamics at i. 
, ) _ E {M /n^(--Qexp(-5) UtS nt+lJlt 
" ' W E { n t} J IL ^ exp ( - S) lit 6n t+1 ,f» 

In the usual manner |2Q, averages involving conjugate fields iu t describe response 
functions, i.e. perturbations of expectation values wrt external fields, so that averages 
involving nothing but conjugate variables correspond to perturbations of a constant and 
will therefore vanish. Moreover, causality implies that G st , which describes the response 
of the average fraction of defaulted companies to at time t to perturbations at time s 
must vanish for s > t. At the saddle point, therefore, we have k t = 0, \m t = 0, Q st = 0, 
iq st = 0, and G st = for s > t, thus iG st = for s < t. 

With these observations we find that the functions $ and ^> appearing in the average 
generating function (jA.lOj) are zero at the saddle point, 

$ = , ^ = , (A.20) 

and the dynamic action S of (jA.14|) simplifies to 

<S = -]rJ2[{l - p)5 st + J 2 q st iu s m t + J2m ( u t - J m t - aJ 2 ^2G st n s \ (A.21) 

^ St t \ s<t / 

With this form of the dynamic action the system-dynamics is described by an ensemble 
independent effective single-node stochastic process of the form 

nt+x = Ut = n t + (1 - n t )Q ^J m t + aJ 2 E G st n s + ^/pr] - $ + 4>\ (A.22) 

the details of which are self-consistently specified by macroscopic properties of the 
system via the saddle point equations, in that each single site process (i) depends on the 
dynamics of the macroscopic fraction of defaulted nodes m t , (ii) the original Gaussian 
white noise is replaced by a coloured Gaussian noise <pt with correlations depending on 

q s t 



f ) = , (0 s t ) = (l-p)5 st + J\ 



st 1 
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and (iii) a memory term appears in the dynamics, if there is some degree of symmetry 
in the interactions, i.e. if a ^ 0. 

The only site-dependence in the averages (■■■)({) appearing the fixed point equations 
(IA.16|) - (jA.18|) comes from the $j dependence in the update rules f it . By the law of 
large numbers, the sums can therefore be evaluated as an average over the -^-distribution 
in the large N limit 

^B-)«)^/d*(*)(-)(^((-)w>, 

i 

in which (. . .)(#) has the same structure as ()A.19|) . except for the fact that the dynamical 
constraints fu of ()A.7j) are replaced by the f&t of ()A.22j) . The saddle point equations 
thus take the form 

m t = ((nt)^))^ , q s t = ((n s n t )^))^ , G st = ((iu s n t }(#))^ 
Appendix A. 3. Simplification of the Single Node Equation 

The single node equation (jA.22|) is complicated by the fact that it is non-Markovian, 
containing a correlation function coupled to the noise term q st and a retarded self- 
interaction G st . This latter term encodes the physics that a firms performance at time 
t is influence by its neighbours, themselves dependent on the firm itself at times s < t, 
via loops in our network of corporate interactions — in particular short loops arising 
through correlated bi-directional interactions. However, as we argued in section El if a 
firm is bankrupt at time s then the performance of partner firms at time t is irrelevant, 
since the firm will still be bankrupt. In the alternative case when the firm is solvent at 
time t, it is clear from the definitions in the dynamics that it must have been solvent 
at time s < t and thus cannot have affected its partner terms at that time. Thus, the 
retarded self-interaction is zero. 

There is a second simplifying feature in ()A.22jl related to the statistics of the 
coloured noise within our system. On multiplying (jA.22|) on both sides by n s with 
s < t and and first averaging over the noise tj> t = (0i, 2 , • • • , <f>t) one finds (n s n t )(#) = 
{n>min(s,t)) (&) & t fixed since if n s = 1, then n t = 1 due to the absorbing nature of the 
defaulted state, whereas if n s = 0, so is the product n s n t , irrespectively of n t . As a 
consequence we have q st = m min ( Si t), and thus 

(0 s t ) = (1 - p)8 st + J 2 m s , s<t, (A.23) 

Having seen that the memory term in the dynamics vanishes, it transpires that only the 
equal-time version of the noise correlation (<f>t<f>t) = 1 — p+ J 2 m t is required to propagate 
the order parameter mt- One needs 

= (( n t+i)m)^ = m + ^((1 - n t )e(J m t + y/pri + <f>t))^ 

In order to evaluate the average (. . .)(#) over the correlated noise in the second term, 
convert the probability density p(<fr t ) = p(4>i, 4>2, ■ ■ ■ , 4>t) into p(n t ,4>t) - - the joint 
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probability density that the node-variable takes value n t (n t = or n t = 1) and the 
noise variable at time t is in an infinitesimal interval around t ; formally one can write 
this as 

p(nt,<k) = j &4>t-lP(<l>t)KtM®A-i) = I d ^*-l p ^*-l'^*) ,5 nt,nt(*,0 t _ 1 ) 

where n t {d, 4> t _i) is the value of n t for the specific $ under consideration, and a given 
noise history 4> t -i- Writing the joint probability in terms of a conditional as 



p(n t ,(f>t) = p{nt\<j>t)p(<p: 



and noting that p(n t \<p t ) must be independent of the conditioning by causality, and 
finally using 

p(n t ) = (n t } { o)5 nul + (l - (n t }(d))Sn t ,o 
for a given •§ one finally obtains 



m t +i 



1 " (nt) 



l + erf 



^2(1- p +Jtmt)) 



(A.24) 



which agrees with the result of our heuristic reasoning in Sec E2 Note that the condition 
c/N — > is not needed in the present argument. 
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